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Abstract. Brenti and Welker have shown that for any simplicial complex X, the face vectors of 
successive barycentric subdivisions of X have roots which converge to fixed values depending only 
on the dimension of X. We improve and generalize this result here. We begin with an alternative 
proof based on geometric intuition. We then prove an interesting symmetry of these roots about 
the real number —2. This symmetry can be seen via a nice algebraic realization of barycentric 
subdivision as a simple map on formal power series in two variables. Finally, we use this algebraic 
machinery with some geometric motivation to generalize the combinatorial statements to arbitrary 
subdivision methods: any subdivision method will exhibit similar limit behavior and symmetry. 
Our techniques allow us to compute explicit formulas for the values of the limit roots in the case 
of barycentric subdivision. 



1. Introduction 

Throughout this paper, we let X be an arbitrary finite simphcial complex of dimension d — 1, 
and we assume that all vectors and matrices will be indexed by rows and columns starting at 0. We 
are interested in roots of the f -polynomial of X, defined as follows. Let denote the number of 
j-dimensional faces of X. We declare that f^^ = 1, where the (— l)-dimensional face is the empty 
face, 0. The face vector, or f -vector, of X is the vector 

/X / nX rX rX \ 

U-li /O 7 • ■ • J Jd-ll- 

Let t denote the column vector of powers of t, {t'^,t'^~^, . . The f -polynomial f^{t) encodes 

the /-vector as a polynomial: 

Much work has been devoted to the study of /-vectors of simplicial complexes, their close relatives, 
the g- and /i-vectors, and the associated polynomials. As it turns out, the entries of these objects 
encode many combinatorial and algebraic aspects of the complex to which they are associated (see 
[UISlI^ for background and further references). 

We focus on a recent result of Brenti and Welker which may initially appear surprising. Let X' 
denote the barycentric subdivision of X, and more generally let X'") denote the n^^ barycentric 
subdivision of X. 

Theorem 1.1. |2J Let X he a {d — 1)- dimensional simplicial complex. As n grows, the d—1 largest 
roots of /^'"' (t) converge to d — 1 negative real numbers which depend only on d, not on X . 

We provide some geometric intuition and motivation for why this result holds. We offer an 
alternate proof of this theorem based on these geometric observations. In the process, we show how 
to compute the d—1 real values for each d. Our first main theorem is; 

Theorem A. Let X be a (d — 1) -dimensional simplicial complex. Then the d—1 largest roots of 
f-^^ \t) converge to d — 1 values which are the roots of a polynomial Pd(t), depending only on d, 
whose coefficients are listed in the last row of the inverse of a particular matrix, Pd . 
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The entries of the matrix Pd, and of its inverse P^^ , are computed in Section [6l Our calculations 
allow us to compute the 'limit roots' thus obtained. In the examples, we observed that these 'limit 
roots' are symmetrically distributed about the point —2, with respect to the Mobius transformation 
X ^ jqrr- Our second main theorem proves this symmetry: 

Theorem B. For any dimension {d~l), the d—1 'limit roots' are invariant under the map x M- jqrj-- 

In fact, more can be said. The existence of a 'limit polynomial' and the symmetry result hold for 
an arbitrary subdivision method, as we show in Theorem 15.51 . 

For barycentric subdivision, this symmetry can be seen through a beautiful algebraic theorem. 
Barycentric subdivision, considered as a map on /-polynomials, induces a function b : Z[t] — > 
as in Sectional We list the values of b on monomials as coefficients in the formal power series in the 
variable x over by defining B : Z[t][[x]] —J' Z[i][[a;]] by B{J2k>Q 9k{t)x'') '■— J2k>o ^{gk{t))x'' ■ 

Theorem C. In Z[t][[x]], barycentric subdivision satisfies the identity 



This paper is organized as follows. In Section[21 we discuss the geometric intuition and motivation 
behind Theorem ll.il In Section|31 we prove Theorem El In Section|31 we prove the symmetry stated 
in Theorem|Bl and prove Theorem[Cj In Section[Sl we extend the symmetry to arbitrary subdivision 
methods. We end with Section |6l where we compute the entries in the inverse matrix P^^ found in 
Theorem [X] as well as all limit polynomials and roots up to the value d = 10. 

Acknowledgements. The authors would like to thank Laura Anderson for bringing Theorem ll.il 
to their attention, Thomas Zaslavsky for interesting discussions, and Dennis Pixton for introducing 
the authors to each other. 

2. Geometric Motivation 

Brenti and Welker's theorem may be surprising at first; there is no dependence on the initial 
complex X, only on the dimension d—1. However, geometrically this makes perfect sense. Barycen- 
trically subdividing a simplicial complex X over and over again causes the resulting complex A^"^ 
to have far more cells than the original A. Because higher-dimensional cells contribute more new 
cells (in every dimension) upon subdividing than lower-dimensional ones, the top-dimensional cells 
begin to dominate in their number of contributions to subdivisions. For example, think of geometric 
realizations so that A^") is a subset of A. Then a randomly chosen cell of A^") should, with higher 
and higher probability as n increases, be contained in the interior of a top-dimensional cell of A, as 
top-dimensional cells contribute far more cells to A^"^ than other cells. 

Each of the f^-i top-dimensional cells of A contributes the same number of cells to A^"^. Since 
these cells eventually dominate contributions from smaller-dimensional cells, the /-polynomial for 
A*^") can be approximated by f^-i times the /-polynomial for the n*'' barycentric subdivision cr^"-* 
of a single top-dimensional cell ad- Since the roots of a polynomial are unaffected by multiplication 
by constants, the roots of /^*"' [t) converge to the roots of /""<* {t) as n increases. 

By definition, the coefficients of f^^ ' (i) records the number of cells of each dimension occurring 
in cr^"'' . The number of cells in each dimension is bounded by a constant times the number of 

top-dimensional cells. Thus, if we normalize /'^'i ' {t) by dividing by the number of top-dimensional 
cells, we have coefficients which, for each fc, record the density of /c-cells relative to the number of 
top-dimensional cells. As this density is positive but strictly decreases upon subdividing, there is a 
limiting value for the coefficient. Thus, there is a limiting polynomial, with well-defined roots. 
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We now formalize this intuition. 



3. /-POLYNOMIALS OF BARYCENTRIC SUBDIVISIONS 

3.1. Barycentric Subdivision and the Matrix A^. To prove Theorem 1X1 we begin by observing 
the effect of barycentric subdivision on /-vectors. One key observation is that barycentric subdivision 
multiplies /-vectors by a fixed matrix, A^, defined as follows. 

Definition 3.1. For i,j > —1, let Xi j denote the number of j-dimensional faces in the interior of the 
first barycentric subdivision of the standard i-dimensional simplex, where by convention A_i^_i = 1 
and Ai^_i is 1 if z = —1 and otherwise. Let A^; denote the (d -I- 1) x (d -|- 1) matrix, with rows and 
columns indexed by the integers —1,0, . . . ,d — 1, whose entry in the i*'' row and j*'* column is Xij: 
Ad := [Aij-]. 



For example, A5 



Lemma 3.2. Barycentrically subdividing a {d — 1)- dimensional simplicial complex X multiplies the 
f -vector by A^; 
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Proof. The faces of X' can be partitioned according to the lowest-dimensional faces of X containing 
them. Each face of X is a simplex of some dimension i, and thus its interior contributes Xij to the 
total number of j-cells of X' (or, if i 
linearity. 

Corollary 3.3. For any n > 0, 



0, exactly one vertex to X'). The claim then follows by 

□ 



/^""=/-A3. 

Thus, to understand barycentric subdivision, we need to understand the matrix F^. We will 
compute the entries in Fd more explicitly in the following two sections, but for now we simply 
observe a formula, also noted in [2 , which follows from Inclusion-Exclusion: 

Lemma 3.4. If j > i then Xi,j = 0. If j < i, then 



fe=0 



k 



□ 



By this lemma, F^ is lower triangular with diagonal entries A^^^ — /°^* = i\. Thus, the eigenvalues 
of are 0!,l!,2!,3!,...,d!. 

3.2. Limit Behavior of the Roots. We now turn to the roots of the /-polynomials f^^ \t)- 
Note that some results in this section were exhibited in or follow from [2j. We present here a 

self-contained exposition in order to present an alternative and more geometrically motivated proof 

of (most of) Brenti and Welker's result, Theorem ll.il 
By Corollary [331 

/^""(t) = /^ASt. 
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As the greatest eigenvalue of is d\, we normalize f'^^"\t) by dividing by (d!)" - let p^{t) denote 
the result: 

Note that this normalization does not alter the roots. It will also often be convenient to reverse the 
order of the coefficients oi p^{t), with the effect of inverting the roots of Pn (t) (that is, the roots of 
f^^ '(t)) about the unit circle in the extended complex plane: 

q^{t) ■.^t'p^it-'). 

We are interested in the behavior of the roots of p^ (t) and {t) as n goes to infinity, so we are 
interested in the powers of A^. To take powers of A^, we diagonalize, 

Ad = PdDdPd\ 

where Dd is the diagonal matrix of eigenvalues 0!, 1 !,..., d! and Pd is the (lower triangular) diago- 
nalizing matrix of eigenvectors. Thus, = PdD^^P^^ . 

Let Dd :— ^Dd- Let t denote the column vector t in reverse order, t = . . .f^)^. For any 

simplicial complex X, we thus have the following equations: 



(t) = f'PdD-.p-h 
■■{t) = {f^Pd){bd)\p-,')t 

'■(t) = {f''P<i)[Dd) {P, 



Pn' 



The goal of Section IH] will be to describe more precisely the matrices Pd and -P^^. As the 
eigenvalues of A^ are 0!, 1 !,..., d!, for large n, DJJ is dominated by its last diagonal entry, (d!)". In 
the limit, the powers of the matrix Dd = ^Dd converge to the matrix 

"0 • • ■ 0' 

Md,d-^ [ ••. : 

• • • 1_ 

Thus, as n grows, the polynomials p^ and approach the polynomials 

P^it) if'' Pd)Md,d{Pd')t and 

q^it) := {f''Pd)Md.,d{Pd')t, 

respectively, in the sense that each sequence converges coefficient-wise in the vector space of poly- 
nomials of degree at most d. 

By CoroUarv 13.31 and Lemma l3^ we know the leading and trailing coefficients of p;^(t) and 9^(i): 

p^{t) = (d!)-"t'' H h fl^^ and q^ (t) = (d!)~" H h fl^^f^. Hence, in the limit, p^{t) does not 

have as a root, but has degree less than d (one root of the p^ diverges to — oo), while q^{t) is of 
degree d with as a root. Because the polynomials qn (t) converge coefficient-wise to the polynomial 
q^{t) of the same degree, their roots also converge: 

Lemma 3.5. 7 Let (P„(i))„ be a sequence of monic polynomials of degree d that converges to a 
monic polynomial Poo(i) of the same degree d. Then the roots of Pn{t) may be numbered as ri , . . . ,r^ 
and the roots of Poo (t) as r^ , . . . , r'^ in such a way that for all j = I, . . . , d the sequence r" converges 
to r°° for n —> CO. 
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Since the roots of (t) converge to the roots of (t) , it follows that the roots of (t) converge 
to the roots of P'^{t) (with one of the roots 'converging' to — oo). 

Because the matrix is lower triangular and M^.d has only one nonzero entry in position 
{d — l,d— 1), we have 

where ed is the c?*'* unit vector, and cx,d is a constant depending on and Pd- As both and Pd 
do not depend on the amount of subdivision n, the roots of and q^ do not depend on the value 
of cx,d, and thus do not depend on any coefficient of f^. This leads us to the following definition. 

Definition 3.6. Define the limit p-polynomial by 

Pdit) e'^p-\ 

and the limit q-polynomial by 

qd{t) := e'^P^'t. 

In this section we have proven: 

Tlieorem 3.7. The following facts hold: 

(1) The roots of f''^' \t) are equal to the roots ofp^(t). 

(2) The roots ofq^(t) converge to the roots of qd{t), and depend only on the dimension of X . 

(3) The roots ofp^(t) converge to the roots ofpd{t), and depend only on the dimension of X . 

(4) The coefficient of t^ in the polynomial pd{t) is the [d — i)*"^ entry in last row of P^_^^. 

(5) The coefficient of t^ in the polynomial qd{t) is the {i — 1)*'' entry in the last row of P^^^. 

This proves Theorem [X] Note the first two facts give an alternative proof of Brenti and Welker's 
result, Theorem II ■![ except for the fact that the roots are all real (that the roots are negative would 
then follow from the fact that all coefficients of these polynomials are positive). 

In Section |6l we will explore the final two facts of Theorem 13. 71 bv computing the entries of P^^ ■ 

4. Symmetry of the roots 

Our goal is now to show that the limits of the roots satisfy the symmetry stated in Theorem IB] 
We will prove this symmetry for the roots of qd instead of pd, as it becomes a mirror symmetry 
instead of a Mobius invariance. 

Theorem 4.1. For every d, 

qd{t) = {-lfqd{-l-t). 
In particular, the roots of qd{t) are (linearly) symmetric with respect to 

To prove this theorem, we start by examining the subdivision of a single closed simplex. The 
following lemma uses the usual difference operator A on a sequence, which takes a sequence {an}n>o 
and returns the sequence {a„ — a„_i}„>i. We abuse notation by using A{an}n>o to also denote 
the first term in this sequence, with context determining whether the result is a single term or a 
sequence. 

Lemma 4.2. Let as be a closed simplex of dimension ,s — 1. The f -vector of the barycentric 
subdivision a'^ of Us is given by 

ff^ = A''-nr=(/)}, = A'^-^"{(i + iy}i. 
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Proof. Lemma 2.1 in [2] states that 



3=0 



j=0 i=Q ^ ' k=0 



Note that the innermost sum is the Stirhng number <S'(i,j) of the second kind (see [B], page 34). 
Reordering this triple summation, we have: 



j=0 k=Q ^ ' 1=0 

= Ei'-^E(-ir'(0^'^'(^) 

s 

= E^'{/'^''(0}^^'^^ 

where in the third equahty we replace k with j ~ k. □ 

Corollary 4.3. Let X be a simplicial complex. The f -polynomial of its barycentric subdivision 
f^ (t) is given by 

d 

/^'w = E^'i/''w}'^'"'- 

The polynomials Pi and are given by 

d d 

{d}:)p^{t)^Y.^'{poi^)ht''' (d!)<zf (t) ^E^'i^o^'a)};*'- 

j=0 k=0 



Proof. These formulas follow easily from Lemma 14.21 by linearity of the difference operator A. □ 

Taking inspiration from the formula for (t) above, we consider barycentric subdivision as a 
function on polynomials in t defined by 

(1) b:Z[t]^Z[tl g(t)K^E^'i5(0};i'- 

fc>0 

(Note that this sum is finite because the iterated finite differences of a polynomial are eventually all 
zero.) 

For a simplicial complex X of dimension (d — I) we have that 

biqfit))^d\qf_,,it). 

The function b is linear, and thus it is given by its values on monomials. It will be convenient to list 
these values as arranged on the 'clothesline' [5] provided by a formal power series in the variable x 
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over the ring Z[t]. We thus consider a function B on the ring Z[i][[a;]] defined as 

k>0 k>0 

Theorem 4.4. (see Theorem\D\) In 'L[t\\\x\] it holds that 

Proof. We expand the right-hand side as a formal power series over x and compare the coefficient 
of ^ therein with the value of = as given in ([T]). We have 



l-(e^-l)i ' ^\'^\m 

^ ' j>0 j>Q ^m=0 ^ 

^E"(t(,l)(-ir"E^ 

j>0 ^™.=0 ^ / fe>0 

= E(E'^E(:)(-ir""'')iJ 

fc>0 ^j>0 m=0 ^ ^ ^ 



□ 



To investigate the stated symmetry, we consider the following map 
(2) i:Z[t]^Z[tl g{t)^g{~l-t). 

Lemma 4.5. The map l is an involution, and it satisfies 

ibi = b. 

Proof. The map t is clearly linear, so it will suffice to prove the claim for monomials. 
It is easy to see that l is an involution. Moreover, ib = bi, as 

6S(e*^) = J = ^ = B(e(-i~*)=") = Bt(e*^). 

^ ' \l-{e^-l)tj 1 - (e^ - 1) ^ J \ J 

The claim follows with term-by-term comparison. □ 

We are now ready to prove Theorem l4.1l 

Proof of Theorem \4-.l\ Barycentric subdivision has the effect on each p- and g-polynomial of mul- 
tiplying on the right by F before the t and t, respectively, and rescaling by dividing by d\. In the 
limit, the limit p- and q-polynomials are invariant under barycentric subdivision up to this scaling, 
so that 

h{qd(t)) = d\qd{t). 

Since the eigenvalues of F are all distinct, qd is characterized by this identity, and by having leading 
coefficient ff_i. 
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Applying Lemma 14.51 we have 

b{q,{-l - t)) = b{L{q,{t))) = i{b{qd{t))) = L{{d\)q,{t)) = d\{q,{-l - t)) , 
and since the lead coefficient of qd{—^ — t) is the claim follows. 

□ 

5. Symmetry for Other Subdivision Methods 

In general, given any polynomial g{t) G we can consider the polynomial tg(i) = g(— 1 — t). 
The coefficient of t'' in g(t) contributes (— times itself to the coefficient of P in ig(t): this 
contribution is up to sign the number of (j — l)-dimensional faces of the {k — l)-dimensional simplex. 
Thus, we can interpret t as a map on formal sums of simplices, as follows. 

We will think of every simplex a G X as a subset of the vertex set of X. Now we can write 

T(la 

We represent the simplicial complex X as the formal sum J^aex °f simplices, each with 

'weight' 1. 

Let us recall some basics about subdivisions of simplicial complexes, pointing to 4 as a reference 
for a more detailed discussion. In the following we will write \X\ for the geometric realization of a 
given simplicial complex X 4, Section 3.1]. 

Definition 5.1 (Compare Section 3.3 of [4j). A subdivision (not necessarily barycentric) of X is a 
simplicial complex X whose vertices are points of |X| and such that: 

(1) For every simplex a oi X there is a simplex a oi X such that a C \cr\. 

(2) The linear map \X\ — !> \X\ mapping each vertex of X to the corresponding point of |X| is a 
homeomorphism. 

We will identify a subdivision of X with the triple {X,X,(j)), where cj) : X ^ X is the function 
associating to each a the smallest simplex a £ X such that a C \a\. 

Now, a subdivision {X,X,(j)) induces a linear map 

b^ : Z[X] ^ Z[Xl CTH^ ^ CT. 

4>{a-)=tT 

In the following definition, we collect together compatible subdivisions in different dimensions, 
calling the result a subdivision method. This is not the most general definition of a subdivision 
method, and our results might hold in greater generality, but we restrict ourselves to the definition 
presented here because in more general subdivision methods can get notationally quite cumbersome 
without adding to the actual idea. 

Definition 5.2. A subdivision method $ is a collection of subdivisions $ := {(cn, cr„, 0„)}„>o 
such that for every map ik '■ <Jk ^ cr„i identifying a fc-face of the standard m-simplex and every 
permutation of m elements tt, the map (jik is the restriction of to TrikicFk)- This ensures that, 
given any simplicial complex X, the complex ^(X), called subdivision of X according to the rule 
$, is uniquely defined by requiring that every n-simplex of X is subdivided as (cr„, (?„, (/>„) G $. 
A subdivision method is nontrivial in dimension n if (f>k is not the identity map for some k < n. 
Clearly if a subdivision is nontrivial in dimension n, then (j)„ is not the identity map. 

Barycentric subdivision is the subdivision method where cr„ — 2'^" and (pnj ~ {j} for every vertex 
j e an- 
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Given a subdivision method in view of the linearity of 6^ for each subdivision, it makes sense 
to write 

cr) = E] bg,(T. 

As with the map b given by barycentric subdivision, for any subdivision method the induced map 
6$ always commutes with the map l: 

Lemma 5.3. For any subdivision method <J>, t6$ — bq,L. 

Before we prove this lemma, we need some properties of the map l. For this paper, the link of a 
simplex tr in a simplicial complex X is the subcomplex consisting of all simplices r in X such that 
cr n r = and cr U r, thought of as subsets of the vertex set of X, is also a simplex of X. Applying 
i to X, we obtain 



^ ^ / ^ ^ ^ -j^^dim cr— dim r j ^ -j^^dimr+1^ 



(tDt 

(3) =E(-l)''"'U(linkr)-l)r, 

TEX 

where x is the Euler characteristic. That the Euler characteristic satisfies this identity can be found 
ingj. 

We now need to characterize how l acts on simplices. We do so by looking at how l affects a 
(not necessarily pure) simplicial homology manifold. If M is an r-dimensional (not necessarily pure) 
simplicial homology manifold with boundary DM and X is a finite simplicial complex such that M 
is (PL-homeomorphic to) the geometric realization of X, then we let [M] denote the formal sum 
of all simplices oi X. In this setting the boundary submanifold dM C M induces a subcomplex 
dX C X, and we take [dM] to be the sum of all simplices in dX. 

Lemma 5.4. For any (not necessarily pure) simplicial homology r -manifold M with boundary dM 
and dimension r > 0, 

m]) = {-^Y^'{m-[dM]). 

Proof. The link of every simplex cr G [M] is of dimension r — dim cr — 1 , and is a homology ball 
or sphere according to whether cr is on the boundary dM or not. If linker is a homology ball, 
x(linkcr) — 1 = 0, and if linker is a homology sphere, x(linkcr) — 1 = (_i)'"-dim(T-i_ Thus, by 
Equation ([3]), 



\r— dimcr — 1 \ ^ i \ ^ / i \dimcT 

— 

ai^dM ^ ^ ere a A/ 

= E {-ir-'<y-i-ir^\m^[dM]). 

agdM 

□ 

Proof of Lemma \5.Sl By linearity, it suffices to prove that /,6$([cr]) = 6$t([cr]) for any simplex cr, 
where [a] = X^tCo- ''' the sum of the simplices contained in the manifold a. Since (f> is a subdivision 
method, 6$([cr]) will also be the sum of the simplices contained in some manifold ^(cr) of dimension 
dim cr. Also, b^{[da]) = [9$(cr)]. 
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The result now follows from Lemma 15.41 

6^,([a]) = 6$((-l)<i--+i(H-[9a])) 

= (-l)d'--+i(6^([a])-5<,([aa])) 

= tbi,{[a]). 

□ 

Theorem 5.5. For any dimension n and any subdivision method $ which is nontrivial in dimension 
n, there exists a unique 'limit polynomial' pn.^{t), such that, for any {d — 1) -dimensional simplicial 
complex X, the roots of converge to the roots ofpn,<s>{t) as k increases. The roots ofpn,<s>{t) 

are invariant under the Mobius transformation x M> . 

Proof. The proof of this theorem is exactly the same as the proof of Theorem 15.51 in the case of 
barycentric subdivision. The key observations there were that (.6$ = 6$i and that there exists a 
unique eigenvector for the maximal eigenvalue of the matrix realizing the effect of subdivision on 
/-vectors. That this eigenvector is unique in general follows from $ being nontrivial in dimension 
n, and is left as an exercise for the reader. □ 

Remark 5.6. Since the above interpretation is on the level of formal sums of simplices, the most 
natural context in which to study it seems to be the Stanley- Reisner ring K[X], defined for any 
simplicial complex X and any field K. A good introduction to these rings can be found in [5] , where 
some properties of the Stanley-Reisner ring of a subdivision of a simplicial complex are explored. 
This brings us to ask the following questions. 

Question 5.7. Is there a (multi-) complex in each dimension whose f -polynomial is related to the 
limit polynomials P'^{t) or q^{t) ? 

More generally: 

Question 5.8. Is there a geometric interpretation of the coefficients or the roots of p^{t) (equiva- 
lently, q^{t))? 

6. Computations 

We finish this paper by computing explicit values for the limit roots up to d = 10. As observed 
in Theorem 13.71 to compute pd we need to compute the matrix P^^- To do so, we first compute a 
more explicit expression for A^. 

Recall that = [Aij]. 

Lemma 6.1. 

(=0 ^ ^ 

Proof. Starting with Lemma [4.21 it follows from Lemma [3.41 that: 
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fe=0 ^ ^ 

fe=0 ^ ^ i=0 ^ ^ 

fe=0 ^ ^ 

r. 

□ 

Now that we have the coefficients of A^, we may compute the coefficients of the diagonalizing 
matrix by computing the eigenvectors of A^. We may then compute using standard inversion 
techniques. Note that it is possible to compute the entries of both Pd and explicitly in terms 
of the entries of A^ and hence in terms of sums and products of integers. To give the reader an idea 
of the numerical consequences of these calculations, we report below the results obtained using a 
standard symbolic computation program. 

For any k > d, the d*'' row of does not depend on k, and gives the coefficients of qd{t)- Thus, 
we present here the matrix P{^^: 




1 





























" 





1 
































1 


1 





























1 

2 


3 
2 


1 


























2 

11 


13 
11 


2 


1 























1 

19 


25 

38 


40 

19 


5 
2 


1 




















132 
10411 


3004 
10411 


45 
29 


95 
29 


3 


1 

















90 


3626 


61607 


245 


385 


7 


1 











34399 


34399 


68798 


82 


82 


2 





15984 


12351860 


7924 


39221 


56 


70 


4 


1 








33846961 


372316571 


18469 


18469 


11 


11 





983304 


119432466 


2296176994 


536193 


919821 


567 


588 


9 


1 





12980789207 


12980789207 


12980789207 


429266 


214633 


71 


71 


2 





1345248918720 


281136722386176 


4358731100 


42780833020 


1335075 


3478503 


1050 


930 


5 


1 


123031432784730871 


123031432784730871 


67808366729 


67808366729 


448471 


448471 


89 


89 
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The roots of pd{t) arc, for d < 10, approximated by: 



d 


= 2 


-1 














d 


= 3 


-2 


-1 












d 


= 4 


-4.1861 


-1.3139 


-1 










d 


= 5 


-8.3642 


-2 


-1.1358 


-1 








d 


= 6 


-16.096 


-1.4706 


-3.1252 


-1.0662 


-1 






d 


= 7 


-30.121 


-4.8761 


-2 


-1.2570 


-1.0343 


-1 




d 


= 8 


-55.208 


-7.5398 


-2.7664 


-1.5661 


-1.1529 


-1.0185 


-1 


d 


= 9 


-99.626 


-11.537 


-3.8404 


-2 


-1.3521 


-1.0949 


-1.0101 


d 


10 


: -177.68 


-17.474 


-5.3206 


-2.5830 


-1.6317 


-1.2315 


-1.0607 



-1 

-1.0057 -1 

We see that these roots are symmetric about the point —2, with respect to the Mobius transfor- 
mation X In other words, if a is a root of Pd,{t), then so is where note —2 is fixed by this 
transformation. The symmetry is more apparent in the (hnear) symmetry about — i exhibited by 
the roots of qd{t), which are the reciprocals of the roots of pd{t). The roots of qd{t) are, for d < 10, 
approximated by: 



d = 


2 


-1 















d = 


3 


-1 


-.5 













d = 


4 


-1 


-.76112 


-.23888 











d = 


5 


-1 


-.88044 


-.5 


-.11956 









d = 


6 


-1 


-.93787 


-.68002 


-.31998 


-.06213 







d = 


7 


-1 


-.96680 


-.79492 


-.5 


-.20508 


-.03320 





d = 


8 


-1 


-.98189 


-.86737 


-.63852 


-.36148 


-.13263 


-.01811 


d = 


9 


-1 


-.98996 


-.91332 


-.73961 


-.5 


-.26039 


-.08668 


d = 


10 


: -1 


-.99437 


-.94277 


-.81205 


-.61285 


-.38715 


-.18795 





-.01004 
-.05723 





-.00563 



Further computations of the coefficients and roots have been carried out by Hiiseyin Ozoguz at the 
University of Bremen. Interesting patterns emerge. For instance, Ozoguz has noted (without proof) 
that the denominators of the reduced fraction representations of the coefficients are all square-free, 
and that certain ratios between roots appear to converge as d increases. We leave the study of these 
phenomena to future work. 

References 

[1] Margaret M. Bayer and Carl W. Lee. Combinatorial aspects of convex polytopes. In Handbook of convex geometry, 

pages 485-534. North-Holland, Amsterdam, 1993. 
[2] Francesco Brenti and Volkmar Welker. /-vectors of barycentric subdivisions. Math. Z., 259(4):849-865, 2008. 
[3] Victor Klee and Peter Kleinschmidt. Convex polytopes and related complexes. In Handbook of combinatorics, Vol. 

1, pages 875-917. Elsevier, Amsterdam, 1995. 
[4] Edwin H. Spanicr. Algebraic topology. Springer- Vcrlag, New York, 1981. 

[5] Richard P. Stanley. Combinatorics and commutative algebra, volume 41 of Progress in Mathematics. Birkhauser, 
Boston, second edition, 1996. 

[6] Richard P. Stanley. Enumerative combinatorics. Vol. 1, volume 49 of Cambridge Studies in Advanced Mathematics. 

Cambridge University Press, Cambridge, 1997. 
[7] Eugene E. Tyrtyshnikov. A brief introduction to numerical analysis. Birkhauser, Boston, 1997. 
[8] Herbert S. Wilf. generatingfunctionology. A K Peters, Wellesley, MA, third edition, 2006. 

[9] Giinter M. Ziegler. Lectures on polytopes, volume 152 of Graduate Texts in Mathematics. Springer- Verlag, New 
York, 1995. 



